% function y = myfunTest(x)
%  y = -100* (x(2)-x(1)^2)^2-(1-x(1))^2;
% end

function [y, dy]=myfunTest(x)
 y  = -100*(x(2,:)-x(1,:).^2).^2-(1-x(1,:)).^2;
 dy = [2*(1-x(1)) + 400*(x(2)-x(1).^2).*x(1); ...
        -200*(x(2)-x(1).^2)];
%  if nargout>1
%  dy = [2*(1-x(1,:)) + 400*(x(2,:)-x(1,:).^2).*x(1,:); ...
%        -200*(x(2,:)-x(1,:).^2)];
 end